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ABSTRACT 

Currently popular models for Type la supernovae (SNe la) fall into two general 
classes. The first comprises explosions of nearly pure carbon/oxygen (C/0) white dwarfs 
at the Chandrasekhar limit which ignite near their centers. The second consists of lower- 
mass C/0 cores which are ignited by the detonation of an accreted surface helium 
layer. Explosions of the latter type produce copious Fe, Co and Ni Ka emission from 
^^Ni and ^^Co decay in the detonated surface layers, emission which is much weaker 
from Chandrasekhar- mass models. The presence of this emission provides a simple and 
unambiguous discriminant between these two models for SNe la. Both mechanisms 
may produce 0.1 — O.GMq of ^^Ni, making them bright 7-ray line emitters. The time 
to maximum brightness of ^^Ni decay lines is distinctly shorter in the M < Mch class 
of model (~ 15 days) than in the Mch model (~ 30 days), making 7-ray line evolution 



^Lawrence Livermore National Laboratory, Livermore CA 94550 USA 

^Department of Astronomy and Astrophysics, University of California, Santa Cruz, Santa Cruz CA 95064 USA 



- 2 - 



another direct test of the explosion mechanism. It should just be possible to detect K- 

shell emission from a sub-Mch explosion from SNe la as far away as the Virgo cluster with 
the XMM Observatory. A 1 to 2 (meter)^ X-ray telescope such as the proposed Con-X 
Observatory could observe Ka emission from M < M^h SNe la in the Virgo cluster, 
providing not just a detection, but high-accuracy flux and kinematic information. 

Subject headings: stars:X-rays - gamma rays:general - radiation mechanisms - stars:supernovae 

1. Introduction 

The luminosity of Type la super novae (SNe la), the brightest stellar explosions, arises in current 
successful models not from the energy of the explosion itself but from the radioactive decay of ^^Ni. 
Following Arnctt (1980), in all current models these explosions are the thermonuclear incineration 
of a white dwarf. The very compact nature of this progenitor ensures that the explosion energy is 
efficiently converted to kinetic energy of expansion (c./. Pinto & Eastman (2000b) for a discussion 
and further references). If a significant luminosity is to be developed, energy must be injected at 
later times when the column depth of the ejecta has declined significantly and radiation can escape. 

It is by now widely recognized that this later deposition of energy results from the radioactive 
decay of ^^Ni. This isotope is a by-product of burning to nuclear statistical equilibrium (NSE), the 
process which is responsible for liberating much of the energy which disrupts the star. While 7-rays 
from ^^Ni decay have yet to be observed from SNe la, the evidence from the temporal behavior of 
the light curve (Colgate &: McKee 1969; Clayton et al. 1969) and from optical and infrared spectra 
(Axelrod 1980; Weaver et al. 1980; Kuchner et al. 1994) is very strong. It is further strengthened 
by the success of the radioactive decay model in explaining the light curve of SN 1987A (Xu 1989; 
Pinto & Woosley 1988b, a; Pinto et al. 1988), from which 7-rays were directly observed (Gehrels 
et al. (1988) and references therein). Just as the X- and 7-ray emission from SN1987A yielded new 
and important information about the dynamics of its explosion, so too observations of SNe la at 
these energies hold the promise of significant advances. 

Despite 30 years of observation and theoretical study, progress towards understanding the 
detailed set of events and the physics leading to SNe la explosions remains elusive. In particular, the 
last decade has seen intense observational efforts in the optical and near IR and matching theoretical 
activity. Though we have learned a great deal about SNe la, it must still be said that incorporating 
current data into a tight theoretical picture remains an enterprise fraught with difficulty. SNe la 
explosion models arc frustrated by uncertainties about thermonuclear flame physics and progenitor 
evolution, while the interpretation of optical data is hindered by uncertainties in the non-LTE 
atomic physics of nearly-neutral silicon-, and especially, iron-group ions. 

If there exists a "standard model" for SNe la, the current definition would be a C-0 WD in 
a close binary system which is pushed by accretion very near the Chandrasekhar mass limit, 1.38 
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M0. As the central density of the dwarf rises above 10^ g cm~^, carbon ignites near the center, 
leading to an outward-propagating thermonuclear flame. This flame is Rayleigh- Taylor unstable 
and very quickly becomes turbulent, leading to a marked acceleration in its progress, as necessary 
to achieve an energetic explosion (c./. Woosley &; Weaver (1986)). This turbulence also makes 
direct numerical simulation intractable (Niemeyer k, Woosley (1997); Khokhlov et al. (1997), and 
references therein); current models are thus hampered by a lack of predictive power. 

A possible alternative to the standard model, or perhaps an addition to it, is a C-0 WD in 
the mass range CGM© — O.9M0, bound to a helium main sequence companion burning helium at 
it's center. Accretion in such systems was studied by Limongi &; Tornambe (1991) who found 
that for accretion rates near ~ 3 x 10~^ Mq — yr~^, of order 0.2 of hehum could be accreted 
until this layer detonated near its base. As others have noted, this is very near the accretion rate 
estimated by Iben k. Tutukov (1991) who found this rate, driven by gravitational wave radiation, 
to be insensitive to the mass ratio of the two components, thereby alleviating the need to fine-tune 
progenitor properties to result in such an explosion. 

The possibility that this ignition mechanism might lead to a Type la supernova was first 
suggested by Livne (1990), and studied in 2-D by Livne & Glasner (1991). In one dimensional 
models, detonation in the helium layer produces an inward moving, focused compression wave 
which drives the central density above 10^ g cm~^ and temperature above 10^ K, causing central 
carbon ignition under explosive conditions. The subsequent evolution and nucleosynthesis was 
studied by Woosley & Weaver (1994) in 1-D, and by Livne & Arnett (1995) in 2-D. While the 
progress of the explosion is rather different in 2-D, both calculations produce very similar results 
which may indicate the robustness of this mechanism. The models possess a number of attractive 
properties, such as production of 0.1 — 0.9 of ^^Ni, with the remainder of the original CO WD 
going to silicon group isotopes. As Woosley & Weaver have pointed out, such sub-Mch models 
may be the production sites for '^'^Ca (produced as ^^Ti), and ^^Ti (made as ^^Cr), neither of which 
are accounted for by either Type II supernovae or Mch SNe la (Timmes et al. 1995). Even more 
attractive is the fact that, since burning takes place at lower densities than in Mch stars, these 
models do not suffer from the problem of excess electron capture and the resultant over-production 
of rare neutron-rich species, such as ^^Fe, ^^Fe, ^^Cr, and ^^Ni which plague Mch models (though 
accretion at extremely high rates accompanied by strong winds may be able to cure this, c.f. 
Brachwitz et al. (2000)). The amount of ^^Ni produced increases monotonically with the mass 
of the WD, implying greater maximum brightness for more massive stars. The increased mass 
may also mean a longer diffusion time, and the correlation of both these properties is in the right 
direction to explain the maximum brightness-decline rate relationship described by Phillips (1993). 
There is some hint of this in the bolometric light curves computed by Woosley & Weaver and by 
Livne & Arnett, but a more accurate comparison with observations requires multigroup transport 
with velocity broadened line opacities as in Pinto &: Eastman (2000b). 

It is not known what fraction, if any, of observed SNe la are due to sub-Mch explosions. Iben 
&: Tutukov (1991) estimated the rate for He symbiotic systems to explode at the rate of 1 per 
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century in the galaxy, which compares well with the SNe la rate determined by van den Bergh 
& Tammann (1991), Cappellaro ct al. (1997), and Hamuy & Pinto (1999). In principle, it ought 
to be possible to distinguish Mch and sub-Mch explosion models from their predicted optical light 
curve and spectral evolution. Years of tuning Mch models to achieve a match between observed and 
computed spectra have resulted in a defining list of properties which the successful la explosion 
must possess. "Successful" models include the Nomoto et al. (1984) Model W7 (Harkness 1991) and 
the Woosley & Weaver (1991) Model DD4 (Kirshner et al. 1993; Eastman 1996; Pinto 1997; Pinto 
& Eastman 2000b). Sub-Mch models have, on the other hand, remained largely unexplored, though 
some sub-Mch models share many of the same desirable properties as their more massive cousins. 
The spectra of current sub-Mch models reported in the literature are too blue (Nugent et al. 1997; 
Hoflich et al. 1996) resulting in part from the presence of iron and radioactivity in the outer layers. 
They also have not produced, to date, Ca at sufficiently high velocities to match a typical SNe la. 
We feel these results are far from conclusive, however. Models for sub-Mch explosions have not 
yet undergone the same degree of tuning as their more massive cousins to bring them into better 
agreement with observation. Current spectrTim and light curve calculations are rendered uncertain 
by their lack of time dependence, especially in the high velocity, low density surface layers. 

Since the sub-Mch models have less mass (e.g. 0.8 versus 1.38 Mq), they might be expected 
to reach maximum light in less time than a Mch explosion. Estimates by Contardo & Leibundgut 
(1998) and Riess et al. (1999) give rise times for nearby SNe la in the range 19 to 23 days in 
B. It would be nice if radiation transport simulations of the light curve evolution were accurate 
enough to conclusively rule out sub-Mch models (or Mch models) by comparison with observations. 
Unfortunately, the effective opacities, and even much of the basic physics, used in most calculations 
performed to date remain highly uncertain. To be confident in estimates of the rise time, one 
requires knowledge of the UV line opacity of low-ionization nickel and cobalt - opacities accurate 
to a factor much smaller than the mass difference between Mch and sub-Mch SNe models. Errors in 
the opacity, which comes principally from velocity-broadened UV lines of nickel, cobalt, and iron, 
translate directly into errors in the rise time and peak brightness (Pinto & Eastman 2000a,b) . Most 
SNe la light curve calculations (e.g. Hoflich et al. (1993); Pinto & Eastman (2000b)) have been 
based solely on lines from the Kurucz list (Kurucz 1991), which likely underestimates the nickel 
and cobalt opacities significantly. We have compared Rosseland mean opacity values computed 
for a pure cobalt composition at p = 10~^^ g cm~^ and T = 20, 000 K using the Kurucz list and 
transition data from the OPAL opacity code (Iglesias et al. 1990, 1992) and find that the OPAL 
opacity value is approximately 5 times larger. Thus, in the Mch calculations by Hoflich et al. (1996), 
the model with the slowest rise time still reaches bolometric maximum light in only 15 days. In a 
light curve calculation by Pinto & Eastman (2000b) of a similar model, the predicted bolometric 
maximum light was at 15 days, even though the B and V light curves did not peak until 20 days. 
However, the available evidence is that the time of bolometric maximum corresponds to the time 
of B maximum. Such discrepancies are enough to raise questions about the accuracy of current 
calculations and to doubt claims that sub-Mch models are, at present, quite ruled out on theoretical 
grounds. 
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In this paper we present a simple test of the sub-M^h model which is quite insensitive to most 

details of either the modeling or the exact nature of the supcrnovac themselves. It is based upon the 
fact that the surface helium detonation in most sub-Mch models produces significant yields of ^^Ni 
at high velocities which are absent in the Mch models. High energy photons produced by the decay 
of ^^Ni can therefore escape largely unimpeded from these surface layers, while photons released 
by decay in the radioactive core are strongly attenuated. A lack of observed high-energy emission 
at early times from SNc la would thus argue strongly, and probably fatally, against the sub-Mch 
model. On the other hand, detection of early emission at high energies would argue strongly in 
favor of these models, as it is difficult to produce significant quantities of radioactivity in the surface 
layers of Mch explosions. (While hydrodynamic mixing in pure C-0 models might conceivably lead 
to significant ^^Ni at high velocities as well, the compositional stratification deduced from early- 
time spectra would be destroyed by such a process.) As we shall show, the decay of surface ^^Ni 
in sub-Mch explosions produces considerable Fe, Co and Ni Ka emission between 6 and 8 keV, 
at flux levels great enough to be detected from extragalactic supernovae, possibly by current and 
upcoming missions, and almost certainly by the proposed Con-X Observatory. They also emit '''^Ni 
decay lines which, due to the 6.1 day half- life of '''^Ni, peak earlier and at higher luminosities than 
in Mch models. 

Most aspects of 7-ray transport in SNe la (and supernovae in general) have been thoroughly 
explored and reported on elsewhere. The 7-ray and hard (Compton scattering) X-ray continuum 
evolution of typical Mch explosion model was discussed by Gehrels et al. (1987), and especially 
by Burrows & The (1990) and Burrows et al. (1991). Clayton Sz The (1991) investigated 7-ray 
transport in Mch SNe la models, specifically W7, and discussed the importance of bremsstrahlung 
emission in forming the keV X-ray continuum. Hoflich et al. (1998) presented results of Monte 
Carlo 7-ray transport calculations for both Mch and sub-Mch explosion models, and pointed out 
differences in ^^Ni 7-ray line light curve evolution (see section 3). However, the present work is 
the first to describe the X-ray properties of sub-Mch explosion models, and to propose observations 
which clearly discriminate between explosion models. In addition, there are significant differences 
between our 7-ray results for sub-Mch models and those obtained by Hoflich et al. (1998) which we 
shall describe below. 

The remainder of this paper is organized as follows: in section 2 we briefly describe the methods 
used for our calculations. In section 3 we summarize properties of the models we have investigated 
- more detailed descriptions can be found in the original references - and present our results, 
describing and comparing the X-ray and 7-ray spectral evolution of Mch and sub-Mch models. In 
section 4 we discuss prospects for positively detecting the unique properties of a sub-Mch event 
with current and future X- and 7-ray missions. In section 5 we describe the evolution of 7-ray line 
profiles and what we might learn with sufficient sensitivity. Section 6 gives a short summary of the 
results. 
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2. Computational Methods 

The transport of nuclear decay 7-rays was computed using an updated version of the Monte 
Carlo (MC) 7-ray transport code FASTGAM (Pinto & Woosley 1987). FASTGAM includes pair 
production and photoelectric opacities and incorporates Compton scattering in the limit of zero 
electron temperature. It follows photons emitted by nuclear decay, pair annihilation, and fluo- 
rescence following K-, L-, and M-shell vacancies (induced both by photoionization and electron 
capture). These photons are followed until they escape from the supernova or are destroyed by 
absorption. FASTGAM computes the energy deposition which results from these processes, the 
primary electron energy spectrum produced by Compton recoil, and the emergent 7- ray spectrum. 

Because the supernova ejecta are ionized, these primary electrons, with kinetic energies rang- 
ing up to ~ 1 Mev, lose energy primarily by exciting collective plasma oscillations. Other loss 
mechanisms include atomic ionization and excitation, and, as first noted by Clayton & The (1991), 
they also copiously produce bremsstrahlung X-rays. Under certain conditions which we discuss 
below, these X-rays dominate the continuum emission from ~ 1 to ~ 50 keV. 

Because the transport equation is linear in the emissivity, the bremsstrahlung spectrum can 
be calculated separately from the MC code and then added linearly to the primary 7-ray spectrum 
calculated by FASTGAM. Therefore, following the MC calculation, the bremsstrahlung contribution 
was computed deterministically using a 1-D, spherical, multi-frequency comoving frame transport 
code (Eastman & Pinto 1993). The opacity is dominated by K- and L-shell photoionization, followed 
by electron scattering. The spectrum of primary Compton electrons, S{E) {e~ s~^ erg~^ atom"^) 
resulting from the MC calculation is used to compute the bremsstrahlung emissivity and solve for 
the emergent X-ray continuum flux. 

The bremsstrahlung emissivity was computed using the continuous slowing down approxima- 
tion, given by 

ri. = ^ r dES{E) V Nz ^ ^^4^(^ dE' [ergs cm'^ s'^ Hz"! ster"!] (1) 

where L{E) is the loss function: the e-folding path length for the electron's energy (c./. Axelrod 
(1980)). The partial bremsstrahlung cross section, daz{E')/dv, for an electron of energy E to 
produce a photon of frequency v while interacting with a nucleus of charge Z, was taken from the 
numerical calculations of Kissel ct al. (1991). The output of the deterministic transport calculation 
is the monochromatic photon energy density as a function of depth through the gas and the emergent 
X-ray flux as measured in the observer frame. 

The X-ray continuum from bremsstrahlung emission is capable of producing additional K-shell 
vacancies. For instance, photons with energy > 7.117 keV are capable of ionizing iron. These addi- 
tional vacancies contribute to the Ka production rate, and were not included in the MC calculation. 
A further source of vacancies comes from direct impact ionization of the K-shell by non-thermal 
electrons, however The et al. (1994) have shown that the collisional ionization contribution is negli- 
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gible at early times. The photoionization contribution to the K-shell vacancy rate was included as 
follows: the monochromatic photon energy density obtained from the deterministic transport solu- 
tion was used to compute the additional K-shell photoionization rate. These were combined with 
fluorescence yields (Kaastra &; Mewe 1993) for Ka and K/3 emission following a K-shell vacancy to 
obtain the line emissivities. These were then used to perform another deterministic transport cal- 
culation (the opacity is the same as for the bremsstrahlung calculation), giving the emergent flux, 
in the observer frame, of K-shell lines produced from photoionization by bremsstrahlung X-rays; it 
is at most a few percent contribution. 

To summarize, the composite spectrum is computed in three steps; 1) the transport of decay 
7-rays is calculated using a Monte-Carlo code which gives the emergent spectrum of unscattered 
and down-scattered 7-rays, the K-shell line flux, and the spectrum of primary, Compton scattered 
electrons. 2) We use the continuous slowing down approximation to compute the bremsstrahlung 
emissivity, which is plugged into a deterministic transport code to obtain the monochromatic X-ray 
photon energy density and the emergent bremsstrahlung x-ray spectrum. 3) Additional K-shell line 
production produced when bremsstrahlung X-rays photoionize Is electrons, is accounted for in a 
final step. All that is needed for this is the photonionization rate, which is computed using the 
monochromatic X-ray photon energy density obtained in step 2. The emergent spectra obtained 
from these three steps are added together to produce the final, composite spectrum. 

The inner-shell photoionization rate due to X-rays is small compared with the rate from pri- 
mary 7-rays. The loss function is determined primarily by the ionization state of the gas, and this 
is determined in turn by the balance of valance-shell photoionization by UV photons and radiative 
recombination. Thus, a fully self-consistent calculation would require a complete solution of the 
radiation transport and statistical equilibrium at all energies from the 7-ray to the infrared. For 
a gas with free (thermal) electron density > nion, however, the loss function is insensitive to 
ionization, a condition which the supernova ejecta satisfy throughout their evolution. Our three- 
step procedure, MC 7-ray transport, bremsstrahlung from primary and secondary electrons, and 
additional K-shell lines from X-ray photoionization, can thus be expected to closely approximate a 
full, self-consistent solution. 

The emergent spectrum computed by FASTGAM consists primarily of velocity broadened 7-ray 
decay lines superimposed upon a smooth, Compton scattering continuum. Sampling the emergent 
intensity with enough energy resolution to produce accurate line fluxes requires a large expenditure 

of computing time. As the opacity for 7-ray lines is absorptive (Compton scattering into the line 
profile, while included in the MC treatment, is negligible), a simpler procedure is to compute 
individual line fiuxes directly using a deterministic transport algorithm. The procedure for this is 
straightforward and was described in Pinto & Woosley (1988a) and Eastman et al. (1994) (we note 
however that that in the present work we have used Ogc = 56° for the mean scattering angle) . This 
deterministic approach to computing 7-ray line transport is approximate, less accurate than MC 
for computing energy deposition, and does not give the Compton continuum; it does do a fine job, 
however, at correctly predicting the i^-integrated emergent line fluxes. This is shown in Figure 1, 
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Fig. 1. — Comparison of the ^^Co 1238 keV line light curve for Model W7, as computed with 
the Monte Carlo code FASTGAM (squares), to that computed using the deterministic transport 
method described in section 2. 

which compares the ^^Co 1238 keV line light curve for Model W7 computed deterministically and 
with FASTGAM. 

3. Results 

The models investigated in this paper consist of two Mch mass models and three sub-Mch 
models. The properties of these models are summarized in Table 1. 

The Mch models include Model W7 by Nomoto et al. (1984) and Model DD4 by Woosley & 
Weaver (1991). The X-ray and 7-ray spectral evolution of Model W7 was studied both by Burrows 
& The (1990) and by Clayton & The (1991). Model DD4 is similar in many respects to Model W7, 
in terms of mass, energy and total mass of ^^Ni produced, although in DD4 the flame speed, which 
in both models was an adjusted parameter, was allowed to exceed the sound speed and become 
a detonation, whereas in Model W7 the flame remained subsonic. The result is that in W7 only 
0.29 Mq of "silicon group" nuclei (by which we mean ^^Mg through ^'^Ca) was produced, whereas in 
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DD4 there is 0.49 M©. Despite this, both models have been shown to give reasonable agreement to 
spectral observations of maximum light supernovac (Harkncss 1991; Kirshner et al. 1993). It is not 
clear whether this indicates a lack a sensitivity of such calculations to distinguish large abundance 
differences, or a large inherent variation in the objects themselves. 

The sub-Mch models investigated include Model 2 by Woosley &: Weaver (1994), and Models 
Ml and M8 from Livne &: Arnett (1995). Whereas Model 2 was a 1-D calculation, Models Ml and 
M8 were 2-D, and have been mapped into a 1-D, angle averaged structure (courtesy of Professor 
Arnett). In terms of the initial C-0 white dwarf mass, the total mass, the amount of accreted 
helium, the total amount of ^^Ni produced and the final explosion energy, these three models span 
a large range of properties. 

Figure 2 shows the composition, density and velocity structure of Model DD4, which may be 
compared to that of Model 2, shown in Figure 3. The biggest difference between the two classes 
of models, for present purposes, is the large mass fraction of ^^Ni on the surface of the sub-Mph 
models. 

Figure 4 compares the 20 day spectra of Models 2 and DD4. One consequence of the surface 
^^Ni in the sub-Mch models is the strong Ka line at 7 kcV. This feature is completely absent from 
Mch models. In Model 2 (and the other sub-Mch models), K-shell line emission is superimposed 
upon a bright bremsstrahlung continuum, dominating the spectrum near 7 keV. At 20 days this 
feature is a combination of contributions from Ni (10 percent) Co (73 percent) and Fe (16 percent), 
reflecting the relative abundances of these ions at the time. The energies and fluorescence yields 
for the 12 contributing lines, taken from Kaastra &: Mewe (1993), are summarized in Table 2. For 
all three elements, roughly 30 percent of the emission is due to the Ka2 transition, 59 percent to 
Kai, 4 percent to KP^ and the remaining 7 percent to K/3i. The effect of the continuum edges of 
the Fe, Co and Ni K-shell photoionization crossections near E > 7 keV is clearly evident in the 
Model 2 spectrum. The jump is strong in Model 2 because the the outer material is rich in Fe 
peak elements. No such jump is present in the spectrum of Model DD4, because the Fe abundance 
is much lower in the outer part of the ejecta, dominated as it is by silicon group elements and 



Table 1: Properties of Explosion Models 



name 


Mtot 


Mco 


Mffe 


E51 


M(56Ni) 


M(56Ni) (surface) 


M(Si group) 


ref 


Model 2 


0.90 


0.70 


0.20 


0.90 


0.43 


0.09 


0.29 


a 


Ml 


0.70 


0.55 


0.15 


0.69 


0.14 


0.03 


0.26 


b 


M8 


1.10 


0.90 


0.20 


0.25 


0.71 


0.17 


0.25 


b 


W7 


1.38 


1.38 





1.2 


0.63 


0.0 


0.29 


c 


DD4 


1.39 


1.39 





1.2 


0.63 


0.0 


0.49 


d 



Table 1: Properties of explosions models considered in this paper. For a more detailed description, 
refer to the original references: a) Woosley & Weaver 1994; b) Livne Sz Arnett 1995; Nomoto, 
Thielemann & Yokoi 1984; d) Woosley & Weaver 1991. 
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M [Msun] 

Fig. 2. — Composition, density and velocity structure of the Mch delayed detonation Model DD4 
(Woosley &; Weaver 1991), at 1000 seconds after explosion. 
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Fig. 3. — Composition, density and velocity structure of tlie sub-Mch, edge-lit helium detonation 
Model 2 (Woosley & Weaver 1994), at 1000 seconds after explosion. 
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unburned carbon and oxygen. 

Although both Models DD4 and W7 have 0.63 M0 of ^^Ni, and therefore produce a large 
number of Ka photons, the optical depth at 7 keV remains large enough to absorb the X-ray 
emission from the ^^Ni core until several hundred days after explosion. This is shown by Figures 5 
and 6, which compare the optical depth at 7 kcV to the "''^Ni mass fraction in Models 2 and DD4, 
respectively, at 100 days. As the opacity is not affected by ionization (for low mean ionization), the 
optical depth at other times can be obtained by scaling as r oc t~'^. Most of the optical depth at 
7 keV is due to L-shell photoionization. Only emission from material at r < 1 is able to escape the 
surface. Consequently, the emergent K-shell line emission is due entirely to the presence of ^^Ni on 
the surface. 

Spectra from Model 2 at various times during the first 160 days are shown in Figure 7. These 
should be compared with the spectra of Model DD4 at the same times after explosion from Figure 8. 
At times t < 100 days, the > 50 kcV spectrum in both cases is dominated by narrow (sec below) 
7-ray decay lines superimposed on a Compton scattering continuum produced from down-scattered 
7-rays. As the column density declines with time, the Compton optical depth decreases below 
unity and the Compton continuum eventually disappears. An excellent discussion of the evolution 
of the Compton continuum is given by Xu (1989) and Xu et al. (1991). The smooth continuum at 
< 50 keV is due to bremsstrahlung emission. 

In Figure 7, the strong feature at 14.4 keV which is visible in the 100 day spectra and afterwards 

is a nuclear decay line of ^^Co. The abundance of ^^Co was assumed to be given by the solar 
^'^Fe/^^Fe ratio (0.027), times the ^^Ni mass fraction in the model. As with the K-shell lines, this 
line will only be visible in a sub-Mch explosion - the optical depth to core ^''Co is too great for any 
appreciable escape to occur. 

The time evolution of the integrated K-shcll line luminosity for all models is shown in Figure 9. 
The effect of the surface ^^Ni is dramatic. Again, because of the large optical depth to core ^^Ni, 
only Ka emission produced on the surface is able to escape, with the result that even Model Ml, 
which has only 22 percent as much ^^Ni as Model W7, is nearly 100 times brighter in the Ka line 
than W7. 

The bremsstrahlung continuum time evolution behaves somewhat differently than that of the 
K-shell lines - at least in the Mch models. Figure 10 compares the time evolution of the integrated 
5 to 8 keV continuum luminosity, excluding contributions from K-shell line emission, for each 
of the models. In the Mch models, the emergent X-ray flux must all come from the core. The 
continuum emission which reaches the surface is produced at somewhat higher energy, where the 
K-shell bound-free opacity is lower. However, the Compton optical depth is not negligible, and 
photons are Doppler shifted by the bulk expansion of the gas as they diffuse out. The amount by 
which a photon's energy is reduced in diffusing to the surface can be estimated as 

AE I m^pK 

— — — ~ (2) 

E ct ct ^ ' 
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Fig. 5. — This figure shows the optical depth at 7 keV (sohd line) and ^^Ni mass fraction (dotted 
line), versus mass, in Model 2 at 100 days after explosion. Only the gas situated at optical depth 
of order unity contributes to the emergent K-shell line flux. 

where I is the total path length traversed, R is the distance between point of emission and the 
surface. Using RK,Vt and p ~ 3M/47^i^■^ this can be expressed as 



10 



M 



02 



15 days\ ' 



8 X 10^ cm - s~^ 



(3) 



E \1AMqJ VO.2/ V * / V ^ 

It is this Dopplcr shifting of photons as they diffuse out which accounts for the fact that there is 
E ^ 7 keV X-ray continuum emission in Mch models. The photons are emitted at higher energy, 
where the K-shell bound-free optical depth is lower, and therefore they have an enhanced chance of 
making it to the surface and escape. Additionally, since the abundances in the outer layers of the 
Mch models are of lower Z, the bound- free optical depth in these layers is less than in the sub- Mch 
case; more X-ray continuum photons from the core will make it to the surface. 

In the sub-Mch models, the presence of a large iron peak abundance in the outer layers ensures 
that the optical depth from K-shell photoabsorption is very large, preventing core X-rays from 
emerging at the surface. The 5-8 keV X-ray continuum in this case is due only to surface ^^Ni. For 
the sub-Mch models, the controlling factor is the K-shell bound-free optical depth in the surface 
layers. Roughly, 

^5-8 keV oc [S56{t) (1 - exp(-r-y))] exp(-rx), (4) 

where -S'56(t) is the 7-ray emission rate ([ergs g~^]), r-y is an effective optical depth for trapping 
7-rays, and Tx is the optical depth the 5-8 keV range. The term in square brackets is a decreasing 
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Fig. 6. — Same as Figure 5, but for Model DD4. In this case the ^^Ni (^^Co and ^^Fe at 100 days) 
sits at large optical depth, and K-shell photons produced in the core are absorbed before they are 
able to reach the surface. 



function of time, but the cxp(— rx) is an increasing function of time and dominates the light curve. 
Thus, for the sub-Mch models, the flux increases with time because more of the surface ^^Ni is 
"exposed" by the declining absorptive optical depth. For the M^h models, we can write 

poo 

L5-8 keV OC [S56{t) (1 - exp(-r^))] exp(-rx') / <j>{E') dE' , (5) 

JE 

where tx' is the Compton optical depth at X-ray energy, in the absence of photoabsorption, (^{E) 
is a normalized X-ray emission distribution function, and 

Er^E{\ + — ^). (6) 

Equation (5) expresses (very approximately) the fact that the surface flux at energy E results from 
emission at higher energy in the interior. At the earliest times, when the Compton optical depth is 
high, the Doppler shift incurred by a photon as it random walks its way out is substantial, and a 
large fraction of the X-rays are absorbed before arriving at the surface. As the optical depth drops, 
there is less Doppler shift and less absorption, and more of the X-rays can escape. Once again, the 
time derivative of the square bracketed term in Equation (5) is negative, while that of the second 
term is positive. 

It is possible that vigorous hydrodynamic mixing in a Mch explosion could bring enough ^^Ni 
up to the surface that the X- and 7-ray light curve would more closely resemble that of a sub- 
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Fig. 7. — Spectral evolution of Model 2. 

Mch explosion than a Mch explosion. As mentioned in the introduction, we regard such extensive 
mixing as unlikely because the compositional stratification deduced from early-time spectra would 
be destroyed. However, for the sake of argument, we consider the most extreme case possible, which 
is complete homogenization of the composition of Model DD4. This is shown in Figure 9 as Model 
MDD4. Unlike either the unmixed Mch models or the sub-M^h models, the 5 to 8 keV light curve 
in this case is nearly flat, and comparable in brightness, at 200 days, to the faintest of the sub-Mch 
models considered, Model Ml. The constancy of the light curve is evidently due to the deposition 
term in eqn (4) decreasing at the same rate that the exponential attenuation term is increasing. 

Below, we discuss prospects for observing the X-ray emission from SN la with current, up- 
coming, and proposed X-ray observatories. In the case of the Chandra Observatory, the sensitiv- 
ity begins to fall off dramatically beyond 5 keV and at 7 keV is quite small (effective aperture 
~ 100 cm^). This precludes direct observation of iron peak K-shell emission from all but very 
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Fig. 8. — Spectral evolution of Model DD4. 

nearby SNe la. We consider the possibility of observing a SN la at lower energy with Chandra; 
Figure 11 shows the 1-5 keV luminosity evolution for each of the models. Here wc find that the 
faintest of the sub-Mch models, Model Ml, is indistinguishable from the two Mch models during 
the first 150 days. However, from 150 days onward, the light curve flattens out and remains nearly 
constant for several hundred days after explosion, whereas for the M^h models it drops off, with an 
e-folding time of ~ 130 days. Off course, model Ml would be very faint optically as well. 

Most previous work on the 7-ray evolution of SNe la focused on the evolution of Mch models 
(Gehrels et al. 1987; Burrows Sz The 1990), and especially on the ^^Co line evolution, since these 
have the greatest chance of detection from an M^h SNe la. However, in the sub-M^h models the 
predicted ^^Ni line fluxes are as much as 10 times brighter than for models such as W7 and DD4, 
within a factor of two as bright as their ^^Co lines. The ^^Ni 158 keV light curve for all models 
is shown in Figure 12. Once again, because the Compton optical depth to core ^^Ni remains high 
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Fig. 9. — Evolution of the integrated K-shell line luminosity ( Ka2, Kai, K/?3, K/3i). 



during the first few °°Ni half-lives, the emergent flux in the sub-Mch models is dominated by decay 
of surface ^^Ni. Not only are the sub-Mch models brighter, but they also peak substantially earlier 
(~ 10 days) than the M^h models (~ 30 days). Similar results are obtained for all the other ^^Ni 
decay lines. Figures 13 and 14 show the light curves for the ^^Ni 750 keV and 812 keV lines, 
respectively. 

Hoflich et al. (1998) have also shown that ^^Ni 7-ray lines peak earlier and are brighter in 
sub-Mch explosions. Indeed, it appears from their Figure 8 that in the two sub-Mch models they 
studied the ^^Ni lines reach maximum luminosity at time t = 0. It is not clear how this can be so. 
In the models studied in this paper, the exploded accretion layer has sufficient Compton optical 
depth to delay the maximum to ~ 10 days. Hoflich et al. (1998) also predict a substantially higher 
luminosity in ^^Ni lines. For instance, their Model HeD6 is halfway between Model 2 and Model M8 
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Fig. 10. — Luminosity evolution in the 5 to 8 keV bremsstrahlung continuum, excluding contribu- 
tions from K-shell line emission. 

in its initial mass but has nearly the same amount of surface ^^Ni (0.08 Mq) as Model 2 (0.09 Mq). 
For the ^^Ni 812 keV line they predict a maximum luminosity of ~ 8.5 x 10"^^ photons s^^. But 
for Model MS, which has 0.17 M© of ^^Ni on the surface, we obtain only ~ 3 x 10^''' photons 
at maximum, 10 days after explosion. We cannot explain this discrepancy; to achieve their line 
fluxes, the surface layers must expand with a much higher velocity than in Models 2, Ml, or M8. 
We note that our line transport results have been verified by performing the calculations by both 
Monte Carlo and deterministic methods, as described in the last section; consistent results were 
obtained between the two methods in all cases. 

The situation is quite different for ^^Co lines. The light curves for ^^Co 847 keV and 1238 keV 
are shown in Figures 15 and 16, respectively. Since ^^Co is so much longer-lived than ^^Ni (ti/2 = 
77.1 days versus 6.1 days), the light curves for these lines reach maximum once the Compton 
optical depth to the core is less than 1, which happens at around 50 days past explosion. The Mch 
and sub-Mch models are distinguishable by the rise time, which is earlier in the sub-Mch models, 
but this might be difficult to detect, particularly for objects which are only marginally within 
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Fig. 11. — Luminosity evolution in the 1 to 5 keV continuum. 

detection limits at maximum brightness. Differences in the core ^^Ni mass can mask the effect of 
any contribution from surface ^^Ni. 

4. Prospects for Detection 

Measurement of the 7 keV K-shell emission from Type la supcrnovae provides a direct and 
straightforward basis for discriminating between the two current classes of progenitor models. De- 
tection of a large K-shell line flux or bright bremsstrahlung continuum would indicate surface ^^Ni, 
which at present is predicted in significant quantities only from the sub-Mch model for Type la su- 
pcrnovae. Additionally, the light curves of ^^Ni 7-ray decay lines from sub-Mph models are distinctly 
different than for Mch models, providing an unambiguous indication of the presence or absence of 
^^Ni at the surface. 

Unfortunately, prospects are at best marginal for detecting X-ray emission from a SN la at a 
distance as great as Virgo using the current generation of X-ray observatories. XMM, which was 
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Fig. 12. — The ^^Ni 158 keV light curve for all models. The sub-Mch models are up to ~ 10 times 
brighter, and peak substantially earlier (~ 10 days) than the Mch models (~ 30 days). 

launched 10 December 1999, has the largest effective aperture at 7 keV of current and near-term 
missions. We have used the XMM simulator (SCISIM) to assess the likelihood of detection for the 
X-ray spectrum of Model 2 at 100 days. Assuming a nominal distance of 15 Mpc, the integrated 
5-10 keV flux is 1.82 x 10"^ photons s'^ cm'^. Using both EPIC MOS CCDs and the PN CCD, 
and assuming a half-power diameter of 15", 10 counts would be detected in a 10^ second exposure. 
Although this may seem like a small number, XMM's background count rate in the 5-10 kcV range is 
very low, 1.7 x 10~^ counts s~^); the estimated background is 1.7 counts, making this (formally) 
a 7a detection. Although this will not provide the the type of detailed kinematic information which, 
ideally, one would like to have to confirm the origin of the counts, it is nonetheless sufficient to 
distinguish with high confidence between the M^h and sub-Mch models for Type la supernovae. 

The prospects for detecting a sub-Mch SNe la with the Chandra X-ray Observatory are lower: 
at 7 keV the effective aperture of Chandra is roughly 100 cm^ for the ACIS-S imaging spectrometer; 
under the previous assumptions this would give only 2 photons in a 10^ second exposure. The 
sensitivity is higher at lower energy, but the supernova is also fainter. The total photon flux 
from Model M8 at 15 Mpc is 3.4 x 10"^ photons s"^ in IkeV < E < 2keV, and 1.5 x 10"^ 
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Fig. 13.— Same as Figure 12, but for the ^^Ni 750 keV line. 

photons s~^ in 2kcV < E < 5keV. The effective aperture of the ASCA Imaging Spectrometer in the 
IkeV < E < 2keV range is approximately 600 cm^. Consequently, a 10^ second observation would 
result in only 2 counts. For the 2keV < E < 5keV range, the effective aperture is approximately 
300 cm^ and the number of counts is approximately 4.5. Both of these count rates would be lost 
in the noise, which is substantially greater. 

The proposed Constellation-X Observatory (Con-X: see http://constellation.gsfc.nasa.gov) of- 
fers the best hope for high resolution spectroscopy of K-shell emission from a sub-Mch SNe la in 
Virgo. The proposal calls for Con-X's effective area to be 15,000 cm^ at 1 keV, and 6,000 cm^ 
at 6.4 kcV. Using the assumptions above, the number of counts received in 10^ seconds would be 
110 photons. Con-X might also be able to detect the lower energy bremsstrahlung continuum. In 
Model 2 at 100 days, the integrated 1-5 keV continuum flux at 15 Mpc would be 2 x 10^^ s^^. 
Taking the mean aperture over this energy range to be 11,000 cm^, the number of counts in a 10^ 
second exposure is 21 photons. 

Immediate prospects for detecting ^^Ni 7-line emission are less promising. As an example, 
the peak 812 keV luminosity for model M8 is ~ 2 x 10^^ photons s"-*^ (Figure 14). At 15 Mpc 
this corresponds to a flux of 7.4 x 10""^ photons cm^ s"-*^, which is just at the sensitivity limit for 
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Fig. 14.— Same as Figure 12, but for the ^^Ni 812 keV line. 

INTEGRAL to detect narrow lines in a 10^ second observation (Winkler 1998). As Timmes & 
Woosley (1997) have pointed out, flux levels such as produced by the present set of sub-Mch models 
could easily be seen with the proposed ATHENA 7-ray Observatory (Johnson et al. 1995). 

The view of Hoflich et al. (1998) for detecting ^^Ni in sub-Mch explosions with INTEGRAL was 
slightly more optimistic than ours. However, as described in the proceeding section, their predicted 
line luminosities are much higher than ours, peaking almost immediately after explosion. 

5. 7-Line Profiles 

While neither current nor planned 7-ray spectroscopy missions are sufficiently sensitive to 
obtain useful profile data from other than a very nearby supernova, line profiles may some day 
serve as another possible probe of the radial distribution of radioactivity. 

Figure (17) shows the evolution of the ^^Ni 1.562 MeV and ^^Co 1.772 MeV line profiles. 
While these are weak lines, they are representative, and their proximity in energy makes them a 
good illustration of the difference in evolution of the ^^Ni and ^^Co emission. Because Compton 
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Fig. 15. — The ^^Co 847 keV light curve for aU models. Although the sub-Mch light curve rises 
more quickly than does the Mch light curves, the differences in maximum brightness depend only 
on the total mass of ^^Ni produced. 

scattering through any appreciable angle removes photons from the profile altogether, the line 
shapes reflect the velocity distribution of radioactive material at low optical depth to the observer. 

In the Mch models, radioactive material is only found deep within the supernova. Initially (15 
days), the radiation which escapes is from material closest to the observer (7-rays emitted from 
material on the far side, away from the observer, see too large an optical depth, and are absorbed). 
The peak intensity of the line is therefore blue-shifted by ~ 6000km s~^, and the line is very weak 
as a consequence of the large optical depth. There is also a broad wing to the red resulting from 
Compton scattering. As the ejecta expand and the column depth decreases, a larger volume of the 
ejecta become visible and the red side of the profile increases in intensity. By 100 days, the ejecta 
are nearly transparent (optical depth to the center at 1 MeV is unity) and the profile is that of an 
expanding optically thin sphere, symmetric about the rest energy of the line. This is clearly seen 
in the ^^Co profile, which reaches its maximum intensity at about this time. The short half-life of 
^^Ni ensures that by the time the profile has shifted back to the rest energy it has disappeared. 

In the sub-Mch models, the high velocity ^^Ni layer on the outside leads to ^^Ni and ^^Co 
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Fig. 16. — Same as Figure 15, but for the Co 1238 keV line. 

lines with rapid rise times (c./. Figures 12 and 15). At the earliest times (15 days), the only 
optically-thin paths are those from the side approaching the observer. The resulting line centers 
are blue-shifted by 7000 km with wings extending up to 20,000 km s"-*- to the blue. As the 
ejecta expand (25 days), the optical depth through to surface layers expanding on the opposite side 
of the supernovae becomes smaller, and the red wings of the profiles increase in intensity. While 
the core radioactivity is becoming visible, it makes only a small contribution at this time. By 50 
days, however, emission from the core dominates the profile, which is now virtually identical to 
the Mch model but for the low intensity but very broad wings from the outer layer. In this model, 
the ^^Ni lines are proportionately much stronger at early times due to the low optical depth to 
the surface layers; its short decay time ensures that the ^^Ni lines will always be observed to have 
higher energies than at rest. 

6. Summary 



In this paper we have presented results of X- and 7-ray transport calculations of Mch and sub- 
Mch explosion models for Type la supernovae. We have shown that the X-ray and 7-ray spectral 
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Fig. 17.— Evolution of the ^^Ni 1.562 MeV and ^^Co 1.772 McV line profiles. The top panel shows 
the evolution of these lines for the sub-Mch model hedtbll at four different times; the bottom panel 
shows results from the Mch model DD4. 



evolution of sub-Mch models is distinctly different than Mch models. The presence of surface Ni 
in sub-Mch supernovae would make them extremely bright emitters of iron peak K-shell emission, 
visible for several hundred days after explosion. K-shell emission from a bright sub-M^h located in 
the neighborhood of Virgo would be just above the limit for detection by the XMM Observatory. 
Detection by CHANDRA is unlikely for a SNe la as distant as Virgo. 

Likewise, the ^^Ni 7-ray light curves all display a substantially different behavior in sub-Mch 
and Mch models. Because of the presence of ^^Ni on the surface, the former are much brighter and 
peak much earlier than in Mch supernovae (~ 10 versus ~ 30 days). The ^^Ni 7-ray line emission 
from a bright sub-Mch explosion at 15 Mpc would be just at the limit for detection by INTEGRAL. 

Moderate resolution spectroscopy of 7-ray lines would allow studying the evolution of line 
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profiles. This would provide the most detailed and readily interpreted data on the distribution of 

radioactive elements and column density in SNe la. Unfortunately, there are at present no plans to 
construct an instrument with sufficient aperture to observe any but the very closest (and rarest) 
SNe la in such detail. 
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Table 2: K-shell Lines 



transition E [keV] yield 



Fe Ka2 


6.3915 


0, 


.1013 


Fe Kai 


6.4047 


0, 


.2026 


Fe K/33 


7.0567 


0, 


.0127 


Fe K/3i 


7.0583 


0, 


.0254 


Co Ka2 


6.9151 


0, 


.1084 


Co Kqi 


6.9295 


0, 


.2168 


Co K/?3 


7.6472 


0, 


.0136 


Co K/?i 


7.6489 


0, 


.0272 


Ni Ka2 


7.4611 


0, 


.1226 


Ni Kai 


7.4782 


0, 


.2451 


Ni K/33 


8.2623 


0, 


.0154 


Ni K/3i 


8.2642 


0, 


.0309 



Table 2: Fe, Co and Ni K-shell lines produced by sub-Mch SNe la (taken from Kaastra 8z Mewe 
1993). 



